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BLACK HOLE - BULGE RELATIONSHIP OF POST-STARBURST QUASARS AT Z ~ 0.3 
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ABSTRACT 

The Mbh — a * relation has been studied extensively for local galaxies, but to date there have been 
scarce few direct measurements of stellar velocity dispersions for systems beyond the local universe. 
We investigate black hole and host galaxy properties of six "post-starburst quasars" at z ~ 0.3. 
Spectra of these objects simultaneously display features from the active nucleus including broad 
emission lines and a host galaxy Balmer absorption series indicative of the post-starburst stellar 
population. These are the first measurements of cr* in such objects, and we significantly increase 
the number of directly-measured non-local objects on the Mbh ~ f* diagram. The "post-starburst 
quasars" of our sample fall on or above the locally defined Mbh — cr* relation, a result that is 
consistent with previous Mbh — cr* studies of samples at z > 0.1. However, they are generally 
consistent with the Mbh~ Lbuige relation. Futhermore, their location on the Faber- Jackson relation 
suggests that some of the bulges may be dynamically peculiar. 

Subject headings: galaxies: active, galaxies: evolution, galaxies: kinematics and dynamics 



1. INTRODUCTION 



The M BH - cr * re l ation dFerrarese fc Merrittl 120001: 
Gebhardt etldl [20001: iGiiltekin et al.l 120091: IW00 et al.l 



20101) has been extensively studied in both active and 
inactive galaxies at low redshifts. The existence of 
this and other black hole - host galaxy scaling r elations 
(jMagorrian et al.lll998l iMcLure fc Dunlop|[2002l ) implies 
that the growth of black holes is intimately connected 
with that of the host galaxies. However, the mechanism 
by which this growth is regulated is still an open matter. 

One propose d mechanism is through black hole accre- 
tion feedback. iDi Matteo et al.1 (|2005l ) show that only a 
modest amount of the accretion energy is required to shut 
down star formation in the host galaxy. Their simula- 
tions of galaxy mergers reproduce the local Mbh — 0* re- 
lation with only ~ 5% of the AGN luminosity being ther- 
modynamically coupled to the surrounding gas. On the 
other hand, the need for accretion feedback is not clear. 
IPengl (2007) argues that mergers between galaxies (ma- 
jor and minor) naturally leads to a tight log-linear Mbh- 
Mbuige relation. To investigate the co-evolution of black 
holes with their hosts, it is necessary to measure the re- 
la tion at higher redshifts . 

IRobertson et al.l (|2006T ) study the theoretical evolution 
in the Mbh — cr* relation by simulating galactic merg- 
ers and including both supernova and black hole feed- 
back. They predict that the slope of the relation remains 
roughly constant out to redshift z — 6, while the scal- 
ing of the relation slightly decreases due to an evolving 
Faber- Jackson relation. The effect of this would be that 
cr, increases for a given Mbh as a function of redshift. 
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Conversely, iCrotonl (|2006T ) predicts that the MsH-Mbuigc 
relation evolves in the opposite sense, i.e., black hole 
mass increases with redshift for a given bulge mass. 

Observational investigations into the evolution of the 
Mbh — cr* relation face difficulties. The sphere of influ- 
ence of the black hole cannot be resolved at high red- 
shifts, thus the methods using stellar or gas dynamics 
(|Kormendv fc R ichstonel 119951) to measure Mbh in the 
local universe are unavailable. Alternatively, one can use 
the virial method to estimate black hole masses, which 
invo kes AGN broad line widths and the Rb lr-Lsioq rela- 
tion (jKaspi et alJl200CHBentz et al.ll2009al) . However, in 
type-1 active galactic nuclei (AGN), the active nucleus 
often outshines the host, drowni ng out stell a r features in 
the spectra. Despite this effect, I Woo et al.l (|2006l 120081) 
studied a sample of Seyfert galaxies and found that black 
holes at z = 0.36 and z = 0.57 were over-massive for 
their given ste llar velocity dispersio ns compared to the 
local relation. iCanalizo et al.l ([2012) avoid the problem 
of an over-powering AGN continuum by studying dust- 
reddened quasars between 0.14 < z < 0.37, which have 
spectra that show both active nucleus and host galaxy 
features. They find a result similar to Woo et al. in 
that reddened quasars fall above the local Mbh — c* re- 
lation. The red quasars of Canalizo et al. also fall above 
the Mbh— Lbuigc relation, supporting the notion that the 
black holes are over-massive compared to their bulges. 

Other recent studies have discussed regimes in which 
the local inactive galaxy relation may not hold. 
iGraham et al.l (|201lD calibrate the Mbh — c* relation 
using elliptical, barred and unbarred galaxies and find 
that the slope of the Mbh — cr* r elation can vary ac - 
cor ding to morpholog y. Conversely, IBeifiori et al.l (|2011[ ) 
and IVika et aLT (|2011D find no correlation between black 
hole mass and Sers ic index, n. At the low- mass end of the 
Mbh — cr* relation, iXiao et ail ()201lD find little difference 
in the relation as defined by barred and unbarred Seyfert 
1 galaxies. Xiao et al. also state that inclination angle 
of a galactic disk components may increase the scatter 
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of the M bh — a* relatio n. This is consistent with the 

results of lBennert et al.l (|2011l) . who find that cr* can be 
biased to larger or smaller values by the disk comp o- 
nent of the galaxy. Recently, iMcConnell et al.l (|2011al |bT) 
showed that brightest cluster galaxies have black holes 
that are significantly more massive than predicted from 
the Mbh — cr* relation given their velocity dispersions. 

In this paper we exhibit the first results of our in- 
vestigation of the Mbh — c* relation using a sample 
of AGN whose host galaxies contain luminous post- 
starburs t stellar populations be tween redshifts 0.2 < 
z < 0.4. IBrotherton et all ([2012, (in preparation)) cata- 
log a subset of quasars with strong B aimer absorption 
series visible in their spectra. These "post-starburst 
quasars" (PSQs) are a convenient case for the study of 
the Mbh — c* relation at higher— z, because spectra of 
these objects simultaneously show the broad emission 
lines distinctive of an active black hole as well as strong 
stellar absorption lines indicative of a post-starburst stel- 
lar population in the host galaxies. The moniker "post- 
starburst quasars" is not a definitive classification of the 
objects' luminosities. While some are luminous enough 
to be considered qua sars, others are not. To be consistent 
with previous work (jCales et al.l 120111 : IBrotherton et all 
2010), we will continue to refer to the sample as "post- 
starburst quasars". 

In this paper we present the stellar velocity disper- 
sions (er*) and black hole masses (Mbh) of six post- 
starburst quasars as measured from spectroscopy per- 
formed with the Low Resolution Imaging Spectrograph 
(LRIS) on Keck I. Each of these PSQs were observed 
with HST ACS F606W as part of a snapshot program 
(10588, PI Brotherton, M.). These images and a detailed 
morphologi c al stu dy of the host galaxies are presented in 
I C ales et al.l ([201 ID . In our analysis we rely on photom- 
etry reported by Cales et al., and we adopt the same 
cosmology that they use: a flat universe with H D = 73 
km s _1 , (Im = 0.27, and = 0.73. The paper is orga- 
nized as follows: we discuss the sample and observations 
in section 2, describe our spectral fitting code and a* 
measurements in section 3, and discuss virial black hole 
mass calculations in section 4. We present our main re- 
sults and discussion of possible biases in sections 5 and 
6, respectively. 

2. SAMPLE AND OBSERVATIONS 

We present results for six post-starburst quasars from 
an overall sample of 609 PSQs s elected from the Sloan 
Digital Sky Survey (SDSS) DR3 (|Abazaiian et al.l[2005h 
that show both strong Balmer series absorptions and 
broad emission lines. The f ull details of th e sam- 
ple selection are described in IBrotherton et al.l 120121 
(in preparation)) , but see also ICales et al.l ([20111 ) arid 
IBrotherton et all |2010). Briefly: objects in the sam- 
ple have total Balmer equivalent widths H tota ; = EWh<5 
+ EW H8 + EW H 9 > 2 A. He was excluded from this 
measurement, because it can be blended with Ca II. The 
Balmer break flux ratio was defined by two 100 A regions 
starting at restframe wavelengths 3740 and 3985 A, and 
the selection criterion was set at f °lf ° > 0.9. 

J 3985A' ^T^A 

The continuum S/N in the SDSS spectra were > 8 as 
measured between 4150 and 4250 A. 
We observed six post-starburst quasars of the sam- 




FlG. 1.— HST ACS F606W images of our sample (Proposal 
ID 10588; PI: Brotherton, M.). North is up and East is to 
the left in each of the images. The LRIS slit was placed at 
parallactic angle with a 1" width. Top: 0030 - 1035, 0057 + 0100, 
and 0237-0101. Bottom: 2102+0005, 2118+0056, and 2306-0100. 

pie using t he L o w Resolution Imagin g Spectrometer 
(jOke et all 119951 : iRockosi et all I2010Q at Keck I in 
September 2010. Conditions were excellent with ~ 0'.'6 
seeing. We employed both the red and blue side cameras 
of LRIS to obtain broad spectral coverage. With this cov- 
erage we are able to observe the Balmer series (except for 
Ha and H/3) of our science targets. We observed several 
individual stars to construct host galaxy templates on 
nights in May 2010 and September 2010 using the same 
instrumental setup as for our science targets. 

We observed all of the sample targets with the same 
setup with LRIS blue and red sides. On the blue side 
we used the 600/4000 grism, and on the red side we used 
the 1200/7500 grating. This setup covers observed frame 
3100 A to 6250 A. Because of the small dispersion in red- 
shift of our science targets, this resulted in slightly differ- 
ent rest-frame wavelength coverage. With this setup we 
observed the quasar continuum and Mg II emission line 
on the blue side and the Balmer series on the red side. 
For one object, SDSS 0237 — 0101, we also observed a sec- 
ond red region out to observed frame 7800 A. For each 
target we placed a 1" observing slit across the nucleus 
of the object and oriented it according to the parallac- 
tic angle. In Fig. [J we show HST ACS F606W images 
of our sample with the LRIS slit ovcrplotted (proposal 
10588; PI: Brotherton, M.). Table □ shows the redshift, 
exposure times, SDSS r mag, the aperture size used to 
extract the spectra, and the morphology of each science 
target. 

We followed standard data reduction procedures using 
the NOAO package in IRAF. First we removed bias lev- 
els using the overscan region and cosmic rays from the 
2D sp ectra using the L. A. Cosmic routine ([van Dokkuml 
1200 If) , and corrected for flat fields. We derived a 
wavelength solution from Hg, Ne, Ar, Cd, and Zn arc 
lamps and shifted individual spectra taken through- 
out the night based on sky lines. We fit sky back- 
ground with a low-order polynomial and subtracted the 
fit from the 2D spec tra. Based on the measurements 
of ICales et al.l ([201 ID . we extracted each spectrum us- 
ing an aperture width equal to twice the effective ra- 
dius of the host's bulge componet (i.e., center ±R c ff). 
This is consistent with what has been done in pre vious 
studies (|Gebhardt et al.ll2000t iBennert et all 120111 ) and 
avoids added uncertainty from aperture corrections when 
measuring stellar velocity dispersions. We used three 
standard stars observed throughout the night at paral- 
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TABLE 1 

Observations 



N 


Object 




redshift 


cxp time (s) 


r [mag] 


ap radius (") 


morphology 


1 


SDSS 003043.59 


- 103517.6 


0.294 


1200 


18.41 


2.52 


bulgc+disk 1 , barred-spiral 


2 


SDSS 005739.19 


+ 010044.9 


0.253 


1200 


17.91 


2.28 


bulge+disk 1 , spiral 


3 


SDSS 023700.31 


- 010130.4 


0.344 


2400 


18.19 


0.22 


bulgc+disk, spiral 


4 


SDSS 210200.42 


+ 000501.9 


0.329 


2400 


18.66 


0.98 


elliptical 


5 


SDSS 211838.12 


+ 005640.6 


0.384 


1200 


18.36 


0.44 


elliptical 


6 


SDSS 230614.18 


- 010024.4 


0.267 


1200 


18.76 


0.26 


bulge+disk 1 



Note. — September 2010 Observations. SDSS 0237-0101 was observed with two setups on the red side of LRIS to obtain broad spectral 
coverage. The blue side was exposed twice for 1200 seconds each for a total of 2400 seconds on the blue side, and 1200 seconds i n each of the 
red si de setups. The sixth column shows the extraction aperture radii in arc seconds. These are the bulge effective radii from Cal es et al.l 
(2011). The last column describes the basic morphology of the object. 
a Objects are in the face-on orientation. 



TABLE 2 
Template Stars 



Star 


Sp Type 


HD 93700 


AO 


HD 94980 


A3 


BD +182402 


A5 


HD 95002 


F0 


BD +182401 


F5 


HD 95958 


F8 


HD 97127 


GO 


HD 95676 


G5 


HD 184850 


G8 


HD 95437 


K0 


HD 230519 


K2 


HD 347822 


K2 


HD 347838 


K5 


HD 96761 


K5 


HD 347857 


K7 


BD +212255 


M0 


GSC 01562 - 00224 


M8 



Note. — In the cases where the same spectral type is repeated 
we used an average of the two stars for the host galaxy templates. 



lactic angle to deriv e a flux calibration with data from 
iMassev et all (|1988l) and lMassev fc Gronwalll (119901 ) . We 
also corr ected for Galactic red dening using the extinction 
values of lSchlegel et al.l (|1998f) calculated at sky positions 
by the NASA Extragalactic Database extinction calcu- 
lator. We used the IR AF task deredden, w hich uses the 
extinction function of iCardelli et al.l (|1989D with the in- 
put extinction values to correct the spectra. 

We recovered the spectral region from 3100 to 7250 
A with the 600 line mm -1 grism and 1200 line mm -1 
grating blazed at 7500 A for each template star. The 
S/N of the template stars spectra was >1000. Table [2] 
lists the stars and their spectral types used to create host 
galaxy templates in this project. We conducted the same 
data reduction procedures on the observed template stars 
as on the post-starburst quasars with the additional step 
of correcting for observed radial velocities. 

The blue and red sides of LRIS have different spectral 
resolutions. The Balmer series fell on the blue side of 
LRIS for the template stars, but was shifted to longer 
wavelengths for the quasars in the sample. Before com- 
bining our spectra and fitting, we convolved the red side 
of the quasar spectra with a Gaussian to match the res- 
olution of the blue side. We achieved an intrinsic reso- 
lution of (7 — 90 km s _1 for the lowest redshift object 
(SDSS 0030 - 1035; z=0.253) at rest-frame 3839 A. We 



achieved resolution of tr= 82 km s 
shift object (SDSS 2118 + 0056; z 
3475 A. 



_1 for the highest red- 
=0.384) at rest-frame 



3. STELLAR VELOCITY DISPERSION 
3.1. Fitting Code 

To measure the stellar velocity dispersion of each host 
galaxy, we follow a direct fitting me thod that has been 
used extensively in previous studies. iBarth et al.l (2003) 
use the direct fitting method to fit the spectr a of BL Lac 
object s. More recently it was employed by Kiao et all 
(120111) to fit a sample of Seyfert 1 host galaxies, and by 
IGreene et al.l (l2010bfl to fit t he spectra of H2O mega- 
masers. iCanalizo et al.l (|2012f) use the method to fit the 
host galaxies of dust-reddened quasars. 

We model the spectra with contributions from both the 
host galaxy and the active nucleus. The host template 
is constructed with a composite of stellar spectra taken 
with the same instrumental setup (see § 3.2 for details) 
and broadened with a Gaussian in order to simulate the 
stellar velocity dispersion of the host galaxy. We include 
a power-law continuum attributed to the intrinsic active 
nucleus spectrum. We also include a Legendre polyno- 
mial that characterizes non-stellar emission in the spec- 
tra, including Fe emission and reddening effects. The 
polynomial makes no assumption regarding the redden- 
ing law [e.g., SMC). Thus, the full model function used 
in the fitting takes the form: 

M{x) = {[T(x) ® G(x)] + C(x)} x P(x) (1) 

where T(x) is the host galaxy stellar template, G(x) is 
a Gaussian, C(x) is a power-law continuum (log(/) = Ci 
+ c 2 *log(A)), and P(x) is a third order Legendre polyno- 
mial. The fitting is performed in velocity space (logA). 
There are a total of nine free parameters: the Gaussian 
velocity dispersion, an offset velocity between the model 
and the data, the power-law continuum scaling and slope, 
the four coefficients of the Legendre polynomial, and the 
relative flux contribution of an "old" and "young" host 
galaxy population. The "young" population template 
represents the post-starburst component separate from 
the older population of the host galaxy (see below). 

The fitting code uses the amoeba routine from Nu- 
merical Recipies (written for IDL by W. Thompson) to 
minimize the \ 2 function. The inputs of the amoeba 
routine are the observed galaxy, the template model to 
fit, and the error for each datum in the observed galaxy 
spectrum. We created artificial errors to weight our data, 
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focusing the fitting region on the Balmer series. To block 
emission lines from the fit, we set the pixel errors to a 
value several orders of magnitude larger than the error 
on the unblocked regions. Thus these emission lines do 
not contribute to % 2 , and are ignored during the mini- 
mization. 

3.2. Templates 

We must be careful in defining a host galaxy tem- 
plate that correctly matches the stellar population of the 
galaxies. A template mis-match from the stellar popula- 
tion of the host galaxy may cause an incorrect measure 
of the velocity dispersion. Mismatch between the host 
stellar population and the template compo sition can be a 
source of error in the dir ect fitting method (|Greene 
l2006HBarth et al.l l20021. Furthermore, the Balmer lines 
of early type stars are intrinsically broadened by the com- 
bined effects of turbulence in the stellar atmosphere, ax- 
ial rotation, and pressure broadening of the energy levels 
within the hydrogen atoms. This intrinsic broadening 
may bias our measured cr* if we were to use a host tem- 
plate that did not include early type stars. As long as 
the templates used in the direct fitting include early type 
stars with luminosity weighted proportions that mimic 
the galaxy stellar population, then the broadening func- 
tion simulates the kinematics of the host galaxy, not the 
intrinsic line widths. 

We set out to create host galaxy templates that 
matched the stellar populations in the host galaxies, in- 
cluding the early- type stars. To do this we collected spec- 
tra of a set of stars ranging across 15 spectral types from 
AO through M8. The spectral range covers 3100 - 7200 

A. 

We simulate single-epoch starburst models by deter- 
mining the luminosity weighted proportions of each spec- 
tral type that contribute to individual models. We pro- 
duced a series of templates to compare to the instanta- 
neous burst models o f Chariot fc B ruzual (private com- 
munication, see also iBruzuall 120071) . The series begins 
with weights favoring early type stars to simulate a young 
stellar population, then progresses by removing the ear- 
liest spectral types and weighting the rest of the stars 
accordingly. When compared with the starburst models 
of Chariot & Bruzual, this simulates a range of starburst 
ages by "evolving" early type stars out of the template. 

We compared the Chariot & Bruzual models to each 
template produced in the series as well as to the indi- 
vidual stars that were included in the templates. Using 
the x 2 statistic as an indicator of the goodness-of-fit, we 
identified nine unique matchings between the Chariot & 
Bruzual models and our templates. We note that we 
only calculated x 2 across the spectral fitting region 3600 
- 4800 A in order to obtain the best match in the re- 
gion of interest, but the templates generally matched the 
Chariot & Bruzual models well at longer wavelengths as 
well. 

Thus, we produced spectral templates that match stel- 
lar populations of ages 321 Myr, 404 Myr, 570 Myr, 640 
Myr, 718 Myr, 904 Myr, 4.0 Gyr, 6.0 Gyr and 9.0 Gyr. 
The 321 Myr model was best matched with a single AO 
star. The 404 Myr model was best matched with a single 
F0 star. The 6 Gyr and 9 Gyr models were best matched 
by G5 and G8 stars, respectively. The remaining models 



were best fit by templates that included all spectral types 
from our library with proportions weighted by luminosity 
toward early type stars. 

3.3. Fitting Procedure 

While the host galaxies show a strong post-starburst 
component, they also contain an older stellar population. 
The observed line widths can be fitted using the Gaus- 
sian function, but the line strengths, ratios, and contin- 
uum shape will be determined by the age of the post- 
starburst population and its relative contribution to the 
overall flux. For this reason, we adopt a procedure that 
tests each of the templates that we created and by doing 
so determine an approximate age of the starburst pop- 
ulation. The best fitting templates were combinations 
of the 9 Gyr template (a single G8 star) with templates 
that were < 1 Gyr old. We allowed the relative contri- 
bution of the "old" and "young" populations to be a free 
parameter. The final values of the parameters can be sen- 
sitive to their initial values. This is because the x 2 value 
can converge on a local minimum within the parameter 
space. To ensure we obtain a robust result that converges 
on the global minimum, we re-fit the spectra after per- 
turbing the input parameters. We present the measured 
velocity dispersions (cr* ) , luminosity- weighted fraction of 
the starburst component, 95% confidence interval of cr*, 
and the age of the younger starburst template in Table 
® 

To determine the 95% confidence interval of our cr* 
measurements, we re-fit the spectra holding cr* constant 
accross the range ±30 km s _1 from the best fit cr* in in- 
tervals of 5 km s _1 . We inspected the resulting x 2 versus 
cr* plot to ensure that the measured cr* was indeed the 
value that minimized % 2 , and we fit a parabola through 
the values in order to define more precisely the confidence 
interval. The 95% confidence interval is de fined as the 
range in cr* that increases the % 2 value by 4 (jPress et alj 
[20071) . 

Figure [5] shows the best fits obtained for our targets. 
The observed spectra are plotted in black. The best fit 
model is overplotted in red, and the residuals are plotted 
in blue. 

We obtained an additional longer wavelength spectrum 
that extends to rest frame 5800 A for the object SDSS 
0237 - 0101. Wc fit the spectral region around Mg lb 
with the expectation that features from the underlying 
older stellar population would be dominant in this re- 
gion. The exact region was from 5020 — 5550 A with 
one emission line at 5198 A blocked from the fit. Fit- 
ting this region can be complicat ed due to [Mg/FeJ cor- 
relation with vel ocity dispersion (jjprgensen et al.l ri999t 
I Wold et al J 120071) . However, we found identical results 
between fits conducted with and without the Mg lb ab- 
sorption blocked from the fit. Despite some effort in 
matching the templates to the host galaxy and using vari- 
ous weighting schemes, we were unable to derive a robust 
measure of cr* from this region. The best fit dispersion 
using this region was 155 (+19, —18) km s _1 , which is 
smaller than the dispersion measured from the Balmer 
region, i.e., 214 (+19, —5) km s" 1 ). However, we were 
able to fit cr* values that ranged from ~ 145 ~ 210 km 
s _1 while still maintaining visually similar fits. In Fig- 
ure[3l we show example fits for this region using cr* values 
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TABLE 3 
Best Fit Parameters 













95% 


younger population 


N 


Object 




a"* 


SB 


low-high 


age 


1 


SDSS 003043.59 - 


103517.6 


73 


99 


44 - 79 


718 Myr 


2 


SDSS 005739.19 + 


010044.9 


134 


33 


126 - 146 


321 Myr 


3 


SDSS 023700.31 - 


010130.4 


214 


79 


209 - 233 


640 Myr 


1 


SDSS 210200.42 + 


000501.9 


162 


86 


154 - 166 


404 Myr 


5 


SDSS 211838.12 + 


005640.6 


134 


63 


118 - 137 


404 Myr 


6 


SDSS 230614.18 - 


010024.4 


54 


86 


34 - 61 


570 Myr 



Note. — The velocity dispersions are shown under the column labeled cr* and have units km s" 1 . Velocity dispersions reported here 
have not been corrected for the aperture. The column labeled SB are the fractional contribution of the star-burst component to the total 
flux shown as percentages. The lower and upper limits of the 95% confidence interval is shown with units km s — 1 . In the last column we 
report the age of the template used to represent the young stellar population. 
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Fig. 2. — Spectral fitting results. The data are plotted in black, the best fit model in red, and the redsidual is plotted in blue. Each 
spectrum is plotted in the rest frame. 
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Fig. 3. — Spectral fitting results for the Mg lb region of SDSS 
0237 — 0101. We have plotted the best fits for three values of tr»: 
145, 175, and 210 km . The data are plotted in black, the best 
fit model in red. 

of 145 km s~\ 175 km s _1 , and 210 km s _1 . As in Fig- 
ure [2j the data are plotted in black, and the models are 
plotted in red. The residuals for the three fits shown are 
similar to one another. These fits indicate that there are 
likely model degeneracies that are affecting the determi- 
nation of a* in this region. 

4. BLACK HOLE MASS 

Black hole masses can be estimated using the virial 
method under the assumption that the broad line region 
(BLR) gas is virialized. To calculate the virial parameter, 
we require measurements of the gas velocity and orbital 
radius. The velocity is, in principle, easily measured from 
broad line widths {e.g., Mg II, H/3, Ha). Furthermore, 
reverberation mapping studies have shown a correlation 
between the radius of the BLR (Rbl r) and the contin- 
uum luminosity of the active nucleus (jKaspi et al.ll2000l ; 
iBentz et~a l. 2009a). Thus, the observables to employ the 
virial method are the line width and continuum luminos- 
ity. 

To convert the virial parameter (VP = RblrAv 2 / G) 
to a black hole mass, one must include an additional 
/-factor, which ac counts for unknow n geometries in the 
broad line region. I Woo et ail (|2010l ) have measured this 
factor by scaling the virial parameter of reverberation 
mapped AGN to match the Mbh — relation, and found 
/ = 5.2. Th is is consistent wit h values found by previ- 
ous authors (|Onken et al.ll2004T ) that are about 1.8 times 
higher than an /-factor assuming an isotropic BLR veloc- 
ity field. 

In the following sections we describe how we measured 
the broad line widths and continuum luminosities used 
in the virial black hole mass calculation. We measure the 
Ha FWHM, and the Mg II AA2796, 2803 FWHM to use 
as BLR gas velocities. We measure the 5100 A AGN con- 
tinuum luminosity to use as a proxy for the BLR radius. 
In the following we will refer to XL\ as L3000 or L5100 
for A = 3000 A and A = 5100 A, respectively. Table H 
displays the relevant measurements. Table [5] shows the 
mass estimates and Eddington ratios. 

4.1. Ha 

The H/3 broad emission line is often used to calculate 
black hole masses in single-epoch spectra of AGN. How- 
ever, in our sample, the absorption-line spectrum of the 



host galaxy can significantly contaminate the H/3 emis- 
sion line. Often the broad H/3 line is completely un- 
observed, so H/3 emission is not suitable for measuring 
Mbh- On the other hand, Ha is much less affected by 
the stellar absorption line than H/3, and broad Ha is al- 
ways observed in the spectra of the targets. 

To calculate bla c k hole masses we adopt the relation 
from lGreene et all (|2010al) : 



M BH = (9.7±0.5)xl0 6 ( 



^5100 



10 44 erg s 



-) 



0.519±0.07 



FWHMhh 
Xl 10 3 km s" lJ 



2.06±0.06 



Mr. 



Green e et al. use a virial coefficient from IQnken et al.l 
(2004) (/= 5.5), which is 1.8 times higher than the as- 
sumption of isotropic rando m motions. Their eq uation 
was originally determined by iGreene fc H o ( 2005) , who 
derive a conversion between Ha and H/3 FWHM. Thus, 
the above equation converts the Ha FWHM to that of 
H /3, and uses an /f actor appropriate for H/3. 

IWoo et ail (|201Cf ) measure the Mbh — o - * relation for 
the reverberatio n mapped sample of the Lick AGN Mon- 
itorin g Project ([Walsh et al.l 120091: IBentz et al.l I2009U 
2010). Woo et al. determine the geometric factor for 
/= 5.2^'g by scaling the determined H/3 virial products 
to the local Mbh _ Q"* relation fo r inactive galaxies de- 
term ined bvlGultekin et all (|2009f ) . Because the relation 
from lGreene et all (|2010af) uses the same virial coefficient 
as found by Woo et al., any differences we find between 
our sample and that of Woo et al. will not be due to 
differences between /-factors. 

To measure the widths of the Ha broad emission line 
we used the task specfit in the STSDAS package of IRAF 
(|KrissHl994l ). Our LRIS spectroscopy does not cover the 
Ha region, so we measure the wid th from the SDSS spec- 
tra (DR7; lAbazaiian et al.ll2009f) . The task specfit min- 
imizes the \ 2 using a simplex method and returns the 
central wavelength, integrated flux, FWHM, and skew of 
each line. We fit the region from 6400 - 6700 A using five 
components: a power law continuum, three Gaussians for 
the narrow line components of [N II] A6548, Ha, and [N 
II] A6583, and an additional Gaussian for the Ha broad 
emission line. We held the relative positions of the [N II] 
lines constant for the fitting procedure, and constrained 
all the Gausians to be symmetric. We fit a narrow line 
model using [O III] 5007 on objects where the line was 
visible and scaled the model to the [N II] and Ha narrow 
lines. The specfit task also returns errors on the fitted 
parameters, which are generally large due to the resolu- 
tion and S/N of the SDSS spectra. We show the fits of 
the Ha region for all objects in Fig. |4] 

The stellar Ha absorption line from the host galaxy 
is small compared to each of these components, and sig- 
nificantly more narrow than the broad emission line of 
the active nuc leus. Using spectr a and host galaxy fits 
for PSQs from lCales et al.l ([2012I . (in preparation)) that 
cover Ha, we determined that the host galaxy absorption 
component does not affect the measurement of the broad 
line FWHM. 

The data of I C ales et al.l (|2012t (in preparation)) cover 
the Ha region for the object SDSS 0057 + 0100, and 
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TABLE 4 

Line Widths and Monochromatic Luminosities 









Ha 




Mg II 


L3000 


L5100 


N 


Object 




FWHM 




FWHM mtr 


3000 A 


5100 A 


1 


SDSS 003043.59 - 


103517.6 


4981 ± 428 


1563 


3066 


4 7+1-7 
' -1.2 


1 ,5+0.6 


2 


SDSS 005739.19 + 


010044.9 


4300 ± 210 


1871 


3464 


S6+ ' 07 


0-3ol°;°* 


3 


SDSS 023700.31 - 


010130.4 


7125 ± 770 


2357 


4855 


0.39+Qi 


0.4118:? 


4 


SDSS 210200.42 + 


000501.9 


6132 ± 315 






0.581°;°? 


n 70+0 O6 
u - ' M -0.06 


5 


SDSS 211838.12 + 005640.6 


6303 ± 501 


2278 


2913 


i-4±8:l 


1 7 +0.3 
i -'-0.5 


6 


SDSS 230614.18 - 


010024.4 


3704 ± 719 


1857 


2468 


i-i±S:2 


1 2+ 12 



Note. — Widths are presented with units km s -1 , and luminosities are presented with units xlO 44 erg s — 1 . The Ha FWHM have been 
corrected for the instrumental resolution of the SDSS using R= 2000. The Mg II widths have been corrected for the instrumental resolution 
of LRIS. Instrumental resolution had a small effect on the widths compared to the errors. Errors on the luminosities were determined by 
varying the fitted slope by ±0.1. 



have higher resolution and S/N (~28 pix 1 on the con- 
tinuum) than the SDSS spectrum (S/N~10 pix -1 on the 
continuum). Following the same procedure, the FWHM 
measured using these data agrees within the errors of the 
FWHM measured from the SDSS spectrum. We adopt 
the width measured with the data from Cales et al., be- 
cause the S/N is an improvement over the SDSS spec- 
trum. 

4.2. Mg II 

As with the Ha emission line, there are sev- 
eral prescriptions for calculating black hole masses 
from t he Mg II AA279 6 , 2803 broad emissio n dou- 
blet dMcLure fe Jarvisl l200l iMcGill et all [200l 

iRafiee fe Halll 120111 : IVestereaard fe Osmerl l2009h Of 

these, the calibration from IVestereaard fe Osmerl (|2009T ) 
is based on the largest number of objects and is 
calibrated with L5100 in addition to other continuum 
luminosities. As described in the next section, we 
consider the 5100 A continuum luminosity to be better 
constrained than the 3000 A luminosity, which is 
typically used in conjunction with the Mg II emission 
to calculate black hole masses. The formulation that we 
adopt from Vestergaard & Osmer is: 



M BH = 10 



,, 96 FWHM Mg II ^ 2 , 



lOOOkms 



f) 2 ( 



L, 



5100 



10 44 ergs" 



M (3) 



According to Vestergaard & Osmer, the uncertainty in 
the zero point of this relation is 0.55 dex. This is the 
dominant factor in the uncertainty on the black hole mass 
calculated from the Mg II emission. 

To measure the Mg II width, we first fitted and 
subtracted an Fe II emission template from the spec- 
trum. Iron emission can contaminate the flux surround- 
ing the Mg II emission and affect the measurement 
of the broad line wid t h. W e used the template from 
IVestereaard fe Wilkesl (|200lP i. who derived the template 
from the observed spectrum of the narrow-line Seyfert 1 
galaxy I Zw 1. 

The template by IVestereaard fe Wilkesl (|200lD only in- 
cludes Fe emission that may contaminate the wings of 
the Mg II profile. To account for Fe emission that may 
contaminat e the c ore of the line, iFine et all (|2008[ ) and 
iKurk et all 1)20071 ) add a constant flux under the Mg II 
core equal to 20% of the average Fe flux at 2930 — 2970 
A. This is based on the theoretical Fe emission models of 



iSigut fe Pradhanl ((2003) . To determine if the theoretical 
flux may have affected our Mg II width measurements, 
we subtracted the additional flux from the fit of SDSS 
0237 — 0101 and re- measured the width of the Mg II 
emission. There was no difference measured between the 
two methods, so we do not believe this small amount of 
theoretical flux has affected the measurements of Mg II 
widths. 

We used a fitting routine similar to the one that we 
used to determine stellar velocity dispersions to scale and 
broaden the Fe template before subtraction. The main 
differences between the way that we modeled the Balmcr 
region and the Mg II region are that we replaced the 
host galaxy stellar template with the Fe template from 
iVestergaard fe Wilkesl (|2001l ) and the fitted continuum 
includes contributions from both the active nucleus and 
host galaxy at these wavelengths. We adopted a power 
law to fit the continuum, as we expect that the host 
galaxy contribution to the flux at these wavelengths is 
comparatively small. We fit a total of five parameters: 
normalization and slope for a power law continuum, and 
the scaling constant, velocity dispersion (ap e ), and ve- 
locity offset for the Fe template. 

We first estimated the continuum from a region rela- 
tively uncontaminated by Fe emission: 3010 - 3026 A. 
The result was used as the continuum initial estimate 
for the main fitting routine. The Fe template does not 
include the Mg II emission line, so we blocked the fea- 
ture from contributing to the % 2 during the minimiza- 
tion. This region varied in size for each object depending 
on the strength of the emission. We allowed ape to vary 
up to 2000 km s -1 , but each fit produced a value smaller 
than 1000 km s _1 . The best fit model was then sub- 
tracted from the spectrum leaving a residual spectrum. 
We show Fe template fits to all objects in Fig. [5] 

Objects SDSS 0057 + 0100 and SDSS 2102 + 0005 were 
peculiar, and were not fit well by the Fe template (first 
and second object on the right side of Fig. [SJ) . The spec- 
trum of SDSS 0057+0100 shows some signs of absorption 
to the blue side of Mg II. However, the noise significantly 
increases at the edge of the spectrum, making the con- 
tinuum level difficult to determine. It is possible that 
SDSS 0057+0100 is a low-ionization broad absorption 
line (LoBAL) quasar, but without a clear determination 
of the continuum level at shorter wavelengths, we hesi- 
tate to classify it as such. In order to fit the Fe template 
to this object, we blocked the wavelength region short- 
ward of and including the Mg II emission (< 2840 A) 
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Fig. 4. — Broad- and narrow-line decompositions of the Ha region. The SDSS DR7 spectra are plotted in the restframe in black. The 
models of the narrow emission lines and continuum are plotted in green, and the Ha broad emission line component is plotted in orange. 
The overall fit is plotted in red. The residuals are plotted in blue. The objects in the left column from top to bottom are 0030 — 1035, 
0237 - 0101, and 2118 + 0056. The objects on the right are 0057 + 0100, 2102 + 0005, and 2306 - 0100. 

from the fit, allowing the overall ht to be determined 
by the longer wavelengths. Then we subtracted the en- 
tire fitted spectrum from the region before measuring 
the Mg II width. The spectrum of SDSS 2102 + 0005 
has better S/N on the blue side of Mg II than that of 
SDSS 0057 + 0100, and the profile of Mg II makes it clear 
that the object is a LoBAL. Because the spectrum shows 
strong absorption over the blue wing of the Mg II line, 
we were unable to reliably measure a Mg II line width. 
We do not report a Mg II width measurement for SDSS 
2102 + 0005, but do include one for SDSS 0057 + 0100. 

4.3. Continuum Luminosities 

The continuum luminosity is used as a proxy for the 
BLR radius in the virial method to measure black hole 
mass. We measure the continuum luminosity from the 
power-law continuum of our model. The power-law is 
multiplied by a Legendre polynomial to simulate redden- 
ing effects and non-stellar emission, and the product of 
the continuum and Legendre polynomial represents the 
reddened nucleus free of host contamination. In order to 
obtain the intrinsic continuum luminosities we follow a 
two step process: First, we calculate a scaling factor be- 
tween the product of the continuum and Legendre poly- 
nomial fitted by our spectral code and the HST F606W 
photometry of the nu cleus, which has b een deconvolved 
from the host galaxy (|Cales et al.ll201l"l ). The scale fac- 
tor effectively flux calibrates our fitted spectrum to the 



HST photometry. Then we apply the scaling factor to 
the fitted underlying power law continuum (without the 
Legendre polynomial). 

We measure the monochromatic luminosities from the 
scaled power law. It is important to note that both 
3000 A and 5100 A fall outside the spectral region 
that we modeled. Thus, the luminosites are extrap- 
olations from our fitted quasar power-law continuum. 
We define uncertainty on the continuum luminosities by 
varying the slope of our intrinsic continuum by ±0.1, 
which is the scatter in the average quasar spectral slope 
(|Vanden Berk et alj|200l . The emission at 3000 A has 
contributions from the active nucleus, the host galaxy, 
and iron. W hile we have a template tha t models the 
iron emission (jVestergaard fc Wilkesl 120011 ). attempting 
to separate the three components would include many 
free parameters and model degeneracies may become a 
significant issue. Furthermore, the photometry that we 
use to scale the fitted active nucleus spectrum is in the 
observed HST F606W filter. Because 5100 A is much 
closer to both our spectral fitting region and the F606W 
photometry, we consider L5100 to be much better con- 
strained than Laooo- We present both measured lumi- 
nosities and their formal uncertainties in Table @] The 
true uncertainty for the L3000 measurement is larger than 
the formal uncertainty that we have reported for the rea- 
sons noted above. 
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Fig. 5. — Fe template fitting in the Mg II region. The LRIS data are plotted in black. The overall fitted Fe template is plotted in red, 
and the continuum is plotted in green. The residual spectrum (data - model) is plotted in blue. The objects are in the same order as Fig. [4] 



5. RESULTS 

In the previous sections we have presented our cr* and 
Mbh measurements of six PSQs. We now examine 
the relationship between the black holes and their host 
bulges through the Mbh — cr* and Mbh~ Lbuige relations. 
We further examine the properties of the host galaxy 
bulges with the Faber- Jackson relation to determine if 
the bulges are dynamically peculiar compared to evolved 
systems in the local universe. By examining these three 
relations, we can determine if any observed offset is pri- 
marily due to over/under- massive BHs or dynamically 
perturbed bulges. 

Using the Ha and Mg II emission line widths, we cal- 
culated masses of the black holes in six post-starburst 
quasars and present the values in Table [5] In five of the 
six objects, we measure larger masses based on the Ha 
broad line than the Mg II broad line. We were unable 
to measure the Mbh based on the Mg II broad line for 
SDSS 2102 + 0005, because the line profile for this object 
is affected by an absorption trough on the blue wing, as 
described above. 

The relations from | Greene et al.1 (|2010al ) and 
iVestergaard fe Osmerl (|2009l ) are each calibrated to 
H/3. In principle they should give consistent results, 
because the two prescriptions are measuring the same 
black hole mass, and both Ha and Mg II should 
be reliable indicators of the BLR gas dynamics and 
M BH - That is, according to iGreene fe Hoi (2005), the 
correlation between the Ha FWHM and the H/3 FWHM 



has a scatter on the o r der o f ~ 0.1 dex. Furthermore, 
IVestergaard fe Osmerl (|2009l ) report their calibration 
produces Mg Il-based masses that are consistent with 
the H/3 and C IV-based masses within ~ 0.1 dex. How- 
ever, the la scatter on the zero point of the Vestergaard 
& Osmer relation is 0.55 dex. If we consider this large 
error associated with the Mg Il-based calculation, then 
our calculated Ha and Mg II masses are consistent. 
The discrepancy could arise if, at least for PSQs, the 
geometric /-factor were different between the emission 
line gas species. This is not unlikely, as magnesium has 
a lower ionization potential than hydrogen, and may 
arise from a different location within the BLR. 

5.1. The Mbh — c* Relation 

We have plotted the objects studied here on the Mbh — 
cr* diagram (Fig. Hp. For comparison we pr esent the rela- 
tion measured by iMcConnell et al.1 (|2011a| ). who update 
the plot with compiled measurements from the litera- 
ture. We also inc lude reverberation mapped AGN from 
IWoo et all (|2010D . The best fit from McConnell et al. is 
plotted as a dot-dashed line, and the fit from Woo et 
al. is plotted as the dashed line. We plot the Ha based 
black hole masses in red triangles. The scatter reported 
by both McConnell et al. and Woo et al. is ~ 0.43 dex, 
and the Ha based Mbh all fall within or above this scat- 
ter with respect to either relation. Three of the black 
hole masses that were based on the Mg II measurement 
are consistent with this scatter. 
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TABLE 5 
Black Hole Masses 



Object 



M BH 
Ha 



M BH 
Mg II 


L(, i/L BtW 
5100 A 


Lbol/LEdd 

[O III] 


1.1 


0.037 


0.019 


0.60 


0.023 


0.076 


1.4 


0.009 


0.057 




0.018 


0.013 


1.0 


0.024 


0.027 


0.60 


0.060 


0.062 



SDSS 003043.59 
SDSS 005739.19 - 
SDSS 023700.31 
SDSS 210200.42 - 
SDSS 211838.12 - 
SDSS 230614.18 



103517.6 
010044.9 
010130.4 
- 000501.9 
005640.6 
010024.4 



3.25 
1.06 
3.49 
3.59 
5.70 
1.56 



+3.0 
1.5 
1.1 
0.5 
+4.0 
1.8 
+3.2 
1.7 
+5.0 
2.6 
+ 1.9 
0.8 



Note. - Units of the black hole masses are xlO 8 Mq and were calculated from the equations of Greene et al. (2010a) and 
Vcstcrgaard & Osmcr (2009). Errors on the Mg II-based black hole masses are dominated by the calibration zero point error (~ 0.55 
dex). Eddington ratios assume bolometric correction factors of 10 for the 5100 A luminosities II Woo fc Urry 2002), and 3500 for the [O 
IIL1A5007 luminosities, respectively. 

The objects SDSS 0030 - 1035 and SDSS 2306 - 0100 
appear significantly offset from the Mbh — o~* relation. 



Both of these galaxies were classified by I C ales et al.l 
(2012) as face-on disks. The face-on disk presents a dy- 
namically cold stellar com ponent with small lin e-of-sight 
velocity dispersion, and iBennert et al.l (|201lD showed 
that this orientation can bias the measured v elocity dis- 
persio n toward smaller values. The sample of IWoo et al.l 
(2006) contains four face-on disks, and the objects do 
have lower cr» values than the remainder of their sample. 
However, the magnitude of the bias for their objects is 
only ^30 — 50 km s _1 (Act* ~ 0.1 dex), whereas the 
two objects here have offsets on the order of ~ 140 km 
s _1 (Act, ~ 0.6 dex). Furthermore, as we show below, 
one of these objects is consistent with the Faber- Jackson 
relation. 

One other obj ect in our sample was classified by 
ICales et al.l (|20l2l) as a face-on disk (SDSS 0057 + 0100), 
but it falls within the scatter of the Mbh — c* rela- 
tion. SDSS 0237 - 0101 also has a bulge+disk mor- 
phology, but is not face-on and falls within the scat- 
ter of the Mbh — c* relation. The last two objects 
(SDSS 2102 + 0005 and SDSS 2118 + 0056) are domi- 
nated by bulges and fall outside the defined la scatter of 
the Mbh — o~* relation. 

In Fig. [7] we show the offset of Mbh (top panel) and 
th e offset of a* (bot tom panel) from the relation defined 
bv lWoo et all (2010) as a function of z. Previous studies 
have foun d an offset f rom t h e loca l Mb h — c* relation, 
includ ing IWoo et al.l (|2006l . I2008D and iCanalizo et"aT1 
(2012). The objects in our sample span a redshift range 
that overlaps with the samples of Canalizo et al. and Woo 
et al.. Furthermore, the trends seen by them appear to 
continue in our objects. As discussed by Canalizo et al., 
at first glance this implies very recent evolution in the 
Mbh — c* relation. 

The average Mbh offset for th e sample of all six quasars 
with respect to the relation bv lWoo et al.1 (|2010T l is 1.25 
dex when determined from the Ha broad line. If we 
exclude the object that does not conform to the Faber- 
Jackson relation (SDSS 2306 — 0100, see below), then 
the average offset in Mbh becomes 1.06 dex. This re- 
sult is significantly higher than the scatter in the local 
Mbh — o~* relatio n, wh ich is ~ 0.43 acc o rding t o both 
IWoo et al.l ([27nf)l ) and IMcConnell et~all The 
offset we measure is slightly larger than the average off- 



sets of red quasars fo und bvlCanalizo et al.l (I2012D and of 
non-local Sy Is from IWoo et al.1 (|2006l 120081 ). Note also 
that the Mg II-based black hole masses are systematically 
lower than the Ha based masses. The average Mbh off- 
set, again excluding SDSS 2306 — 0100, based on the Mg 
II measurements is only 0.64 dex. The small number of 
objects in our sample preclude us from claiming whether 
the entire sample of PSQs is on average offset from the 
local Mbh ~ c* relation. To determine that, we will need 
to study more objects from the parent sample of PSQs. 
Furthermore, it is not unreasonable that for these six ob- 
jects we will find some that fall above the ler scatter of 
the local relation, but it is curious that all of the offsets 
are on the same side of the relation. 

The issue is complicated by the matter of the discrep- 
ancy between black hole masses measured with the Mg 
II and with the Ha emission lines. We used the SDSS 
data to measure the Ha emission for all but one object, 
and the SDSS spect roscopy has low S/N (< 10 pix -1 ). 
iDennev et al.l (|2009t ) show that using a Gauss-Hermite 
polynomial to fit low S/N data may overestimate the 
FWHM of the H/3 line. In the one object for which we 
have high S/N LRIS spectroscopy (SDSS 0057 + 0100) of 
the Ha region, we measured a lower FWHM of the Ha 
broad line compared to the SDSS spectrum. The differ- 
ence was ~ 400 km s _1 (~ 8%). The black hole mass for 
SDSS 0057+0100 using the larger width is only a factor of 
1.2 greater than using the smaller width, which is within 
our defined errors. Smaller Ha widths would bring the 
Ha-based mass measurements closer to the Mg II-based 
measurements. However, this small effect is unlikely to 
reconcile the two measurements. As an alternative ex- 
planation, it is possible there is a systematic difference in 
BLR gas dynamics between the Ha and Mg II-emitting 
gas. 

5.2. The Mbh~- -^buigc Relation 

We plot the Mbh— Lbuigo relation in Fig. [HI This serves 
as a test of the bulge properties of the host. If the objects 
follow this relation, then the black holes of these galaxies 
may be "normal" for their bulges, while the velocity dis- 
persions cause the offset seen in the Mbh — er» relation. 
In order to compare the bulge luminosities that we have 
measured here to measurements of systems in the local 
universe, we need to account for the change in luminosity 
associated with the evolution of stellar populations. 
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Fig. 6. — Plot of our objects compared to the locally-defined Mbh — &* relation. The relation by McConncll ct al. (2011a) is plotted 
as a dot-d ashed line and gre y plus symbols. The scatter is shown as a shaded grey region. The relation of reverberation mapp ed active 
galaxies by Woo ct al. (2010) is plotted as a dashed line with grey asterisks. We plot the red quasars from RTanalizo ct al. (2013) in green 
squares. The black hole masses of the PSQs based on Ha and Mg II are plotted as red triangles and blue diamonds, respectively. 
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Fig . 7. — The offset from the local Mbh — o* relation ( W oo et aLl 
2010) is plotted against the redshift of the objects. The PSQs 
from our sample are plotted twice with triangles representing the 
Ha-based Mbh measurement and the diamonds t he Mg Il-based 
measu rement (same as Fig. [BJ. The red quasars of Canali ze- et al.l 
(120121) are plotte d as squares, and th e Seyfert 1 galaxies from 
IWoo et al.l pOOM ) and I Woo et~aTI (|2008l ) are plo tted as X symbols 
The objects that define the local relation from Woo ct al. (2010) 
are plotted as plus symbols. 

We measured the bulge luminosities in a procedure 
that is similar to the measurement of the AGN contin- 
uum luminosities. For each object, we scaled our fitted 
host galaxy temp l ate to the measured bulge photome- 
try of [Cafes et atl pOlTI ) in the HST F606W filter. We 
then measure V band photometry from the scaled host 
spectrum. Photometry from Cales et al. has already been 



k-corrected, but we also apply a p assive evolution correc- 
tion u sing the online calculator^ of lvan Dokkum fc Franxl 
(l200l . We applied the "Default 1: single stellar popu- 
lation" setting using the calculator, which uses a burst 
of star formation that passively evolves. This setting ap- 
proximates the starburst event of our PSQs, which have 
luminosities dominated by the young population (see Ta- 
ble©. 

The "Default 1" setting of the passive evolution cal- 
culator assumes that the stars were formed at redshift 
z » 1. However, the PSQs clearly have had a recent 
starburst within the past < 1 Gyr. Thus, we adjusted 
this parameter to more closely match the star formation 
history of our targets. The difference in redshift from 0.3 
(the average of the sample) to 0.4 corresponds to roughly 
840 Myr. This is appropriate for all of our targets except 
SDSS 2118 + 0056, which has a redshift closer to 0.4 than 
0.3. For this target we adopt the passive evolution cor- 
rection assuming a formation redshift of 0.5 (a 740 Myr 
stellar population). As expected, this produced a larger 
correction than assuming all stars formed at z » 1, 
and applying this correction shifts the objects to smaller 
V band luminosities. We adopt an error on the bulge 
luminosities that assumes 0.1 mag error from the fits of 
Cales et al. and an additional 10% error on the spectral 
host- AGN decomposition. 

In Fig. [8l the solid line shows the best fit relation from 
iMcConnell et al] ()2011af ). who compile measurements of 

4 www. astro . yale .edu/dokkum/evocalc/ 



12 



individual galaxies from the literature, including bright- 
est cluster galaxies. The two dotted lines show fits de- 
fined by the uncertainty in their scaling and slope, while 
the grey shaded region shows the intrinsic scatter they 
measured. Two of the Ha-based black hole mass mea- 
surements (triangles) fall above the scatter in the re- 
lation, while the remaining objects are consistent with 
the relation. However, if we consider the Mg II-based 
black hole masses (diamonds), all of the objects fall on 
the relation. We also note that the scatter from Mc- 
Connell et al. (0.5 de x) is somewhat greater than that of 
iGiiltekin et~aLl (|2009f ). who found 0.38 dex. If we were to 
consider the fit from Gultekin et al., four of the Ha-based 
black hole masses would fall above the relation. Thus, 
to some degree, the interpretation of the offset from the 
relation depends on the relation's true intrinsic scatter. 

The result is quite similar to that from the Mbh 
h , relation. That is: PSQs fall within the scatter or 
above the Mbh— Lbuigc relation if we consider the black 
hole masses measured from Ha. Whereas, if we consider 
the Mbh measured from Mg II, the PSQs are consistent 
with the relation. From this we can surmise that if the 
offset is to be believed, the black holes really are over- 
massive for their bulges. 

5.3. The Faber-Jackson Relation 

Figure [9] shows the location of our samp le on the Faber- 
Jackson relation (jFaber fc Jacksonl [l976'). If the er* val- 
ues that we measured are causing the offset seen from the 
Mbh — c* relation, we expect the objects to fall below 
this relation. That is, the objects should have velocity 
dispersions that are too small for their bulge luminosi- 
ties. This is the case for one object, but the remaining 
five are generally consistent with the relation. 

The left and right panels show the SDSS g band and r 
band absolute AB magnitude, respectively. The g and r 
photometry were measured with the same method that 
we used to measure the V band luminosity (see previ- 
ous section) with the additional step of converting lumi- 
nosities into SDSS AB magnitudes for ease of compari- 
son. This also required transforming the passive evolu- 
tion corrections from B and V to q and r corr e ctions , 
for which we used th e relat ions from lJester et al.l (|2005| ). 
iNigoche-Netro et al.l (|2010[) calculate the Faber-Jackson 
relation in discrete magnitude bins of width 1 magnitude. 
Considering the g absolute magnitudes (left panel), the 
PSQs appear to follow the relations defined by Nigoche- 
Netro et al.. The dashed line represents the relation ap- 
propriate for objects with magnitudes between —19.5 and 
—20.5. Two objects fall in the magnitude bin —18.5 to 
— 19.5 and are consistent with that relation. The only 
object that falls significantly offset is SDSS 2306-0100, 
which is one of the face-on disk galaxies. Considering the 
r magnitudes (right panel), all objects fall in the magni- 
tude bin —20 > M r > —21. Five objects are consistent 
with the relation, while SDSS 2306-0100 is again signifi- 
cantly offset. In both panels we have plotted the overall 
relation using the entire sample of Nigoche-Netro et al.. 
The fit for their overall sample has an intrinsic dispersion 
o f 0.6 dex. 

iDesroches et al.l (|2007l ) have also studied the luminos- 
ity dependence of the Faber-Jackson relation. They de- 
rive a quadratic fit to to the relation in the SDSS r band 
using data from the NYU Value Added Galaxy Cat- 



alog (iBlanton et al.l [2005h . While INigoche-Netro et~all 
(|2010() use only galaxies that are well fit by a de 
Vaucouleurs profile an d the corresponding magnitudes, 
IDesroches et al.1 (|2007t ) select objects with a range of 
Sersic indices and calculate a "Sersic-like" magnitude 
from a Petrosian magnitude. As Desroches et al. explain, 
the Sersic index is correlated with the galaxy luminosity 
and therefore may bias the slope of the relation if not 
accounted for. We have plotted the quadratic fit from 
Desroches et al. log cr* = -1.79+0.674(log L)-0.0234(log 
L) 2 ; M= -2.5(log L)) with 0.1 dex scatter in the right 
panel of Fig. [9] (s olid line and dark gr e y sha ding) along 
with the fit from INigoche-Netro et all (|2010l ). Three of 
our objects fall on the relation derived by Desroches et 
al.. One object falls above the relation (within 2a), while 
two fall significantly below the relation. The two objects 
that are significantly offset are also the most offset from 
the Mbh — relation, suggesting that at least part of 
the offset is due to peculiar dynamics of the stellar pop- 
u lation. 

iNorton et al.l (|200fl) studied the velocity dispersions of 
quiescent E+A galaxies. To separate the post-starburst 
population from the underlying older population of stars 
they simultaneously fit the 4100 - 5250 A region with 
K, G, and A star templates. They find that the post- 
starburst population does not show a significant cor- 
relation between its velocity dispersions and the over- 
all host luminosities. However, the distribution of dis- 
persions for their sample is con sistent with the F aber- 
Jackson relation using dat a from lFaber et al.l (|1989l ) and 
Uorgensen &: Franxl (|1994l ). The older stellar population 
fi* does show correlation with host luminosity, but the 
objects fall on average ~0.6 mag brighter than the ellip- 
tical galaxies of similar dispersion. We do not separate 
the post-starburst population from the older population 
with unique measures of ex* as Norton et al. do, so the 
value of a* that we measure has influence from both pop- 
ulations. 

6. DISCUSSION 

We have measured Mbh and tr* in six AGN that are 
hosted by galaxies showing post-starburst stellar popula- 
tions ("post-starburst quasars"). These objects fall on or 
above the locally defined Mbh — and Mbh— Lbuigc re- 
lations. The location of the PSQs on the Mbh — c* dia- 
gram is consistent with other studies of objec ts at z > 0.1 
(jCanalizo eFaIll2012t IWoo et al.ll200ll2006f ). 

We have used the Balmer region, which includes the 
Balmer series, Ca lines, the G band, and other stellar fea- 
tures to estimate the stellar velocity dispersions. While 
we have used host galaxy templates that approximate the 
stellar population of the host galaxies, it is possible that 
the young and old populations have different velocity dis- 
persions. The kinematics of the young population are in- 
fluenced by the mecha nism that tr i ggered the starburst. 
In a theoretical work, iBekki et~a l. (2005) describe vari- 
ous models of galaxy interactions (major/minor mergers, 
strong tidal interactions, etc.) that lead to E+A spec- 
tral signatures. They find a range of morphologies, ve- 
locity dispersions, and starburst mass fractions resulting 
from the merger event depending on the specific model 
used. The kinematic properties of the resulting post- 
starburst population depend on the initial parameters of 
the galaxy-galaxy interaction that induces the starburst. 
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Fig. 8. — The Mbh - Lb u i go relation. The solid line is the relation from McConncll ct al. (2011a). The dotted lines show the relations 
defined by the uncertainty in the fitted scale and slope with scatter shown as the grey shaded region (±0.5 dex). Ho based black hole 
masses of the PSQs are plotted as triangles, and the Mg II based masses of the same objects are plotted as diamonds. 




Fig. 9. — The Faber-Ja ckson relation. Left: The SDSS g band absolute AB magnitudes. The dashed line is the best fit from 
Nigochc-Nctro ct al. (2010) for the magnitu de bin -19.5 > g > -20.5. Right: The SDSS r band absolute AB magnitudes. The dashed line is 
the best fit from Nieochc-Nctro ct al. (2010) for the magnit ude bin -2 > r > -21. The solid line and dark grey region show the quadratic fit 
and scatter to the Faber-Jackson relation fromjDesrochcs ct al. (2007). For both panels the the light grey shaded region shows the scatter in 
the fit from[Nigochc-Nctro ct aL (±0.25 dex, l2010l ). and the dot-dashed line shows the best fit for the overall sample of Nigoche-Netro et al.. 



The morphologies presented by iCales et al.l ()2011[ ) for 
our sample indicate that these objects have only minorly 
disturbed features, if any, that could indicate a recent 
merger event. Four objects of our sample are described 
by Cales et al. as containing a disk component. However, 
the overall sample of PSQs display a range of morpholo- 
gies, and Cales et al. found that half of the observed 
sample appear to be post-merger remnants. The vari- 
ety of morphologies seen in the parent sample is con- 



sistent with the simulations conducted by iBekki et al.1 
(|2005[ K indicating a range of possible galaxy-galaxy in- 
teractions that trigger the starburst event {e.g., tidal 
interactions, minor/major mergers). Regardless of the 
triggering mechanism, the PSQs, at least in this sub- 
sample, have Mbh and c* that are generally consistent 
with previously studied z > 0.1 galaxies. 

The origin of the velocity dispersion (pressure sup- 
ported versus rotational) should influence whether the 
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objects fall on the Mbh — a* relation, which is best 
denned for pressure supported bulges. Observationally, 
iNorton et all (|200lD studied the young and old popula- 
tions independently in E+A galaxies and found them to 
be largely pressure supported, with the younger post- 
starburst populatio n having higher a* than the old pop- 
ulation. However, iBekki et al.l (|2005f ) had trouble re- 
producing t h is res ult with their models. By contrast, 
iPracv et al.l (|2009D found rotation to play a dominant 
role in their studied E+As. 

We do not separate the young and old populations in 
our analysis, so the kinematical differences between the 
post-starburst population and older population in PSQs 
remains a question. In principle, we could adopt an ad- 
ditional free parameter in our model to measure the ve- 
locity dispersion of the old population separately from 
the young population. However, this would introduce 
degeneracies that would be difficult to disentangle, espe- 
cially with the presence of the active nucleus. We are 
currently conducting an investigation into the kinematic 
differences between young and old populations in quies- 
cent E+A galax ies that lack the active nu cleus continuum 
contamination (jHiner &: Can alizo 20 Q in preparation). 
This analysis will be similar to that of Norton et al., but 
we will include multiple wavelength regions. By includ- 
ing longer wavelength regions, such as the Mg lb and 
Ca II triplet regions, we will be able to further distin- 
guish between the young and old stellar populations. 

We have also measured the bulge luminosities of the 
host galaxies. After plotting the Mbh~ Lbuigc and Faber- 
Jackson relations, we find a scenario which rules out 
cr* bias as the source of the offset from the Mbh — re- 
lation for three of the six objects we measured. Two 
objects (SDSS 0030 - 1035 and SDSS 2306 - 0100) are 
significantly offset from the Mbh — c» relation. SDSS 
2306 — 0100 is additionally offset from the Faber-Jackson 
relation, indicating that the measured <r* may be signif- 
icantly influenced by the face-on disk of the host galaxy. 
SDSS 00 30 — 1035 is consistent with the relation as de- 
fined by iNigoche-Net r o et al.l (120101). but it falls signifi- 
cantly below that of iDesroches et all (|2007[ ). so it may 
also be influenced by the face-on disk. While SDSS 
0237—0101 falls above the Desroches et al. Faber-Jackson 
relation by less than 2a, it is also the object that falls 
closest to the best fit Mbh — &* relation. A reduc- 
tion of ~ 50 km s _1 in ct* would bring it in line with 
the Desroches et al. fit, but drive the object left of the 
Mbh — c* relation (consistent with the other offsets). 
Peculiar dynamics can clearly influence the measured ve- 
locity dispersion, and with a more statistically significant 
sample, the general trends will become more apparent. 

The offset from t he Mbh — o* r e lation is consistent 
with that found by iCanalizo et al. | (|2012f). Simil ar to 
their samp l e and the samples of I Woo etall ([2001 and 
IWoo et all (|2008f) , the black holes measured here all have 
masses > 10 8 Mq. This brings up t he question o f sam - 
ple selection at high-z (z > 0.1). lLauer et al.l ()2007f ) 
describe a potential bias that selects for high-mass black 
holes at higher-z, which is a product of a steep luminos- 
ity function and cosmic scatter in both t he Mb h — c* and 
Mbh - Lbuigc relations. ICanalizo et al.l (|2012f ) also show 
that highly luminous AGN (log(L 5 i o/erg s" 1 ) > 43.6) 
are offset from the Mbh — c* relation regardless of their 



redshifts. All but one of the PSQs in our sample have 
luminosities in this range. 

We have presente d six objects of a l arger class of post- 
starburst quasars (Brother ton et al.1 120121 in prepara- 
tion). This is clearly not statistically significant and we 
may have by chance chosen objects that all fall above 
the Mbh — o* relation. To discuss PSQs as a whole, 
we require measurements of more objects that define a 
statistically significant sample. We have obtained data 
for additional PSQs from the sam e parent sample th at 
will be exhibited in a future work i 1 liner et al.l 120121 in 
preparation). In addition, it will be important to disen- 
tangle the kinematics of the younger post-starburst pop- 
ulation from the underlying older stellar population. For 
this purpose, we have obtained spectroscopy of quies- 
cent E+A galaxies. We will probe possible differences 
in er* kinematics between the populations using both the 
Balmer series and longer wavelength regions (e.g., Mg lb, 
CaT). 
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